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PREFACE 


As indicated in the title, this paper considers the problem of 
the Submarine Commanding Officer who must determine salvo size 
and spacing in making a torpedo attack on a target vessel, 

The paper presupposes a knowledge of certain concepts of 
the theory of probability, The terminology and notation used in 
discussing mathematical probability correspond to the expressions 
in general usage in the literature, such as, ''The Theory of 
Probability", M. E. Munroe, 1951. The paper is offered parti- 
cularly as a consideration for possible classroom discussion by 
opemations analysis students. 

The writer wishes to thank Professor Thomas E. Oberbeck 
of the U. S. Naval Postgraduate School for his original suggest- 
ions, assistance, encouragement and cooperation in the prepara- 


tion of this paper. 
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GLOSSARY OF NOTATION 


= track angle 
= target true bearing 


target true course 


ваш > 
" 


domain or area of uncertainty associated with position of MOT 
f,g = probability density functions 

h = hit probability 

i = aim point designator; i = 1,2,---,n 

L = target length 

a.p.8. = aim point sequence 

n = number of aim points in an aim point sequence 

m = number of hits required for a specified degree of damage 


= number of T/O's fired at a single aim point. 


minimum acceptable probability of a hit 


= target range 


= target speed 


time on station in days 


maximum salvo size 


number of T/O's on board on jth day of patrol; j = 0,1,---,T 


= number of T/O's with which S/M begins patrol 


= expected frequence of attack on jth day of patrol 


e a 2 د وح‎ v dw з 
1 


angle on the bo:v (AOB) 


- maximum expenditure of T/O's in one attack on jth day of patrol 


DM 





CHAPTER 1 
ELEMENTS OF THE PROBLEM 
1.1 Statement of the Problem 

Each time a submarine commander decides to make an attack 
he must make a complete firing phase decision. This consists, 
among other things, of deciding: 

(1) 

(a) How many hits are wanted on the target. 

(b) The price in torpedoes he is willing to pay for these hits. 

Contained also in the firing phase decision is: 

(2) 

(a) Selection of the final inputs into the problem. 

(b) Selection of a torpedo spread technique to insure hitting 

the target. 

(c) Control of the firing of the torpedoes [particularly as 

to time). 

This paper will try to give aid to (1) by an analyais of (2) 
based upon statistical inferences concerning 2(a) and the mathema- 
tics of probability concerning 2(b) and 2(c). 

The basis for decision will be derived in terms of probability 
of hits obtained for torpedoes expended. 

The handling of one attack is then to be correlated with the con- 
duct of a complete patrol based upon what is desired to be acootgiplished 


by this patrol, as specified in the patrol orders. 





1.2 Patrol Objectives 

A patrol is considered as taking station in a geographical area 
under one of two alternatives as to length of stay: 

1. No restrictions as to time, 

2. Stay on station a specified number of days. 

Each alternative is pursuant to patrol orders. The first envi- 
sages pursuing every opportunity to attack, and early expenditure 
of torpedoes relieves the submarine of further commitments in the 
area. The second reflects a need to keep the submarine in the area 
for a definite additional reason, (e.g. reconnaissance, rescue, or 
to provide a constant menace to the enemy over a specified period 
of time, with availability of relief, say, a governing factor). 

It will be considered that submarine patrol orders contain the 
order to "Inflict Maximum Damage on the Enemy"! which may be 
interpreted to specify two clearly distinct tasks: 

A. "Insure" kill of each target attacked, 
or 

B. "Insure" damage to each target attacked. 

Here the word "insure" is used in a probabilistic sense, as are also 
the words "kill" and "damage". 

It is desired that the submarine commander do his best with his 
opportunities in any one of four categories in which his orders place 


him, namely Al, A2, Bl, B2, as shown in the following table. 


l. See Theodore Roscoe, Submarine Operations in World Маг 11, 1949; 
Sub Pac Operational Plan, p. 


$ 


| 











Kill each target attacked. Damage each target attacked. 
Leave when torpedoes Leawe when torpedoes 
expended. expended. 

Space the kills. Space the damages. 

Stay on station a specified Stay on station a specified 
length of time. length of time. 


In a type Z patrol, i.e., when A2 or BZ is assigned, note that it is 
desired that the commander spreadfi his effort over the time speci- 
fied so that his ship presents, in effect, a continuous threat to the 
enemy. [This reasoning should also leave him with some torpedo 
capacity to protect himself in a tight spot throughout the duration 
of a type 2 patrol. ] 

It will be seen that in a type 2 patrol the submarine Command- 
ing Officer must have definite knowledge of the expected frequency 
of enemy target traffic through his patrol area upon which to base 
his estimate of expected frequency of attack. This implies that there 
will be sources of intelligence from which he can estimate this target 
frequency and also that he must have some measure of his own search 
ability. It will be assumed that he can, from previous sources of 
intelligence of whatever origin, (e.g. intelligence reports, previous 
contacts, previous attacks, etc.) make such an estimate of the ex- 


pected frequency of attack on any given day of his patrol. 


275 
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1,3 Firing Phase Objective 


For uniformity, we define the firing phase as the period from 
the time of the last observation (last sef of target data inputs) until 
the time when the last torpedo of the salvo crosses the target track, 
real or extended. 

In any attack, the final firing set-up, determined by whatever 
means, is put into the torpedo data computer at a given instant, 
The input data to be scrutinized here are: True bearing, range, 
angle on the bow(or true course) and speed of the target. True 
refers to the compass rose only; the input data are to be consider- 
ed as a set of random variables in the probability sense. The 
possible errors of these input data are evaluated to predict the 
size and probability distribution of the location error of the target, 
i.e. the nature of an area (or domain) of uncertainty, say D , 
which is moving and expanding with elapsed time, and into which 
the torpedo spread is fired. In particular, it is desired to know 
how unit mass, representing the totality of probability of target 
position, described by the probability density function f(D) asso- 
ciated with D , is distributed over this area D at a certain 
elapsed time t after the final inputs are made. 

From knowledge of D and f(D) at any given elapsed time t, 
it is desired to predict the aim point sequence, spread spacing, and 
salvo size necessary to obtain; (1) at least one hit with a minimum 
acceptable probability Ро when damage is wanted: or, (2) at least 


m hits, with a minimum acceptable probability Бы ‚when m hits, 





on the average, are needed to sink the target. 

The Commanding Officer then completes the firing phase deci- 
sion on the basis of the information cited above. If a type 2 patrol 
is being conducted, it is also necessary to specify the salvo size for 
a given attack that is compatible with the requirements of the length 
of stay on station. The following chapters discuss the elements of 
this decision, i.e. salvo size as a function of the choice of one of 


the four patrol objectives and "how good the set-up is". 





CHAPTER г 
ELEMENTS OF ONE-DIMENSIONAL CONTINUOUS 
PROBABILITY DISTRIBUTIONS 
2.1 Scope 
Certain fundamentals of probability will be discussed, using 
simple models, in order to develop some concepts, definitions 


and notation for later use. 
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272 Target Position 


Consider a target T of length L бп course Со . For сопуеп- 
ience, Со will be taken as the positive direction on the axis of a 
one-dimensional coordinate system. Suppose that, at the instant t, 
the MOT (middle of target) is known to be on the x-axis, but its 
exact position is uncertain, The position of the MOT , at the in- 
stant t , may be regarded as a stochastic variable with an asso- 
ciated probability density function f(x, st). The function f(x, :t) 
will be termed the target position density function (at time t ). 

In the examples considered below, it will be convenient to distin- 
guish between the cases Kite f(x, ;t) is uniform and where 


f(x, +t) is normal. 


Hence, when f(x, :t) is uniform, set f(x, 5t) £ (x,:t) where 


1 1 1 
Ex, for as T3 Ex 5 х За, +5Ех, 
(2.1) f (x, :t) = 
0 elsewhere 


and a, is the value of ху denoting the mean of the distribution. 


(See Fig. 2a). In this context, target location error is taken to be Ex). 


W hen f(x, it) is normal, set f(x, it) = f(x st) where 


1 1 XP To ^ 
(2. 2) f (x,;t) = —— Exp - = 
n 1 У?тбх, 2 | 6x, 


Where -oo< хі < го 
Again, а, is the mean of the distribution but in this context the 


| 
target location error is taken to be 26x, : 


— Torpedo Position 


Let x, represent a point on the same x-axis where a hypothe- 


2 
tical torpedo (hereafter called T/O) crosses the axis at the instant 
t. Then h(x,;t) ‚ the probability of a hit on the target may be 
computed for any point х); i.e. 

X, + 5 L 
(2.3) h(x,;t) = f(x, 5t) dx, : 

x, = = 

J 2 2 

The function h(x, ;t) is shown in Fig. 2 b for the case where the 
target position error Ex, exceeds L. In this connection, set 
Ex, = M*L where M isa constant, М > 1, 

Suppose it is known that a T/O aimed at any point on the x-axis, 
зау х, will definitely cross the x-axis at the instant t , but the 
exact point of crossing is known only uncertainly as a stochastic vari- 
able with the associated conditional probability density function 
g(x\it |х). 

Then, Р(х) , the conditional probability, of a hit on the target 


when the T/O is aimed at the point x4 is given by the formula 


]. Some attributes of the normal distribution are discussed in 
Appendix A 





/Г оо 
(2.44 Р(хұи) = | h(x, it) g(x,it|x,) dx, . 
– со 
By way of illustrating Р(х) » assume glx,;t|x,) is also a 


uniform probability density function such that: 


l 1 
= for хз - ; Ex, < х, 5 ха + z Ex, 
552 
g(x, ;t |x,) = 
0 elsewhere 


~. 


In this context, the quantity Ex, will be termed the T/O error, 
Assume also that f(x, st) is given by (2.1), Finally, let 


Ex > L,i.e., Ex, = MeL , Ex, = NeL , with 


l 2 l Д 
М > М > 1. Тһе graph of P(x, it) under the assumptions of 


> Ex 


this paragraph is shown in Fig. 2 c. 
Whenever the instant t is not essential to the discussion, it 
\ 


will be convenient to write 
f(x, :t) E f(x) 


h(x 


h(x, it) 2) 


P(x, it) P(x 


3) 
eic. 
Since the primary interest of this cbapter is the probability 


density function associated with target position, it will be assumed 


that the T/O error, Ex, » is zero, Under this assumption, 


(2.5) f х, = х, 


|205) - мо) 
2,4 Computation of hit probabilities. 
An hypothetical example is inserted at this point, primarily 


as a vehicle to explain the notation used in the paper, 





o 
AA EE m e c «3 | 
п БЕС 
һы 
x Е ——— 





It is convenient to formulate first of alla system for designating 
the crossing points of T/O's on the k-axis. By assumption [see (2. 5)] 
the crossing point of the T/O is the aim point of the T/O. The choice 
of aim points, and the probable consequences thereof is a basic 
objective of this paper. A definite choice of equally spaced aim 
points will be called an aim point sequence (a.p.8.). Inasmuch as 
it will be assumed throughout this paper that all T/O's cross the 
axis at the — instant t , the a.p.s. is not concerned with the 
order of firing; it specifies only (Бела of aim points, where 
aimed, and the spacing between aim points, Thus, let an a, p. s. 
be defined as follows: 

a, = point marking the mean of f(x) 

d 1/2 distance between adjacent aim points along the axis 

п =number of aim points in the a. P. 8. 

See Fig. 2d. | 


Then, according as n is odd or even the aim points inana,p.8, are: 


n odd n even 
a, = а a, = asta 
а, Fea, + 2d a, = ab td 
а, = а, + 2а az = a, СЕ за 
а, = „АА + 3d 


etc, 
Attention is tur ned fórst to the case where a single T/O is 
placed at each aim point of an a. p. 8., with spacing 2d = L, 


a so-called simple salvo. Later, the case of placing more thẹn one 








T /O at each aim point ofan a.p.s. will be discussed. 


For the case of a simple salvo, let P. denote the probability 


of a hit associated with the ith aim point of ana.p.s., i=1,... ,n. 
n 

Let P (n) = Ж Pi , 
i=] 


1.е. , P(n) is the total probability of one hit for a simple salvo 
of n aim points. 

By way of illustration, Table J.gives values of Р, апа P | (n) 
for simple salvos in two cases. The first case, designated in the 
table by f, , isfor f(x,) uniform with Ex, - 4L [cf. formula 
(2. 1)] and the second case, designated by f, , is for f(x, ) 


normal with Ox, = L [cf. formula (2.2)] 


Table I 
nzl n=2 | n=3 п = 4: 
а.р. fa f fa f fa fa f f 
а, „290.385 2250. 356. 2250 .383 „250 .341 
а, ‚250 Жар 1250 .242 „250 „341 
а, 2-0: 242 73250 29196 
а, 250 „136 


Р „250 „383 .500 „682 .750 .867. 1.000 .954 

In the case of simple salvos, with 2d = L, note that 
there are no gaps nor overlaps in the spread spacing. This 
insures for any n, a maximum value of Р, in the case of sym- 
metrical target position density functions, such as £ x) and 


f x). 


10 








Two possible criteria for the choice of n are as follows: 

(a) Let Ро = minimum acceptable probability of a hit with 
one T/O. In this case, n = maximum i for which Р.Р, . 

(b) Let Po minimum acceptable total probability of one hit 
for a.p.s. selected, In this case, n = minimum i for which 
P,2 P,- 

By way of illustration, suppose Ро is fixed at 0, 2; then from 
Table 1, for f. Say, it can be seen that, n = 4 would not be 
used, but instead n = 3 would be used, If, however, a high 
РЬ, was desired, say 0.9, then n = 4 would necessarily be 
used, The Aa (a) — (b) isa matter of preference. 
Hereaíter, the choice of n will be given in terms of both 
criteria. | 

The possibility of more than one hit with a salvo to T /O's 
arises if 

(1) more than one T /O is fired at each aim point of an 

a.p.s. with the distance between adjacent aim points 
exactly equal to target length, i.e., 2d = Ц; 

or (2) one T/O is placed at each aim point of an a. p. в. 
with the distance between adjacent aim points less 
than target length, iie., 2d < L; 

or (3) more than one T/O is fired at each aim point and 
2d is less than L. 

These possibilities will be considered in the order just enumer- 


ated and will be referred to as Cases 1, 2, and 3, respectively. 
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These cases will be considered under the simplifying assumptions 
that 

(1) each T/O crosses the x-axis exactly at its intended aim point. 

(II) all T/O's of a salvo may be considered to cross the x-axis 

at virtually the same instant t , even if more than one 
T /O is fired at any aim point, 

In considering Cases 1, 2, and 3, let Вп, тп, 2d) denote the 
probability of obtaining two hits with a salvo of size nxm in which 
m T/O's are placed at each of the n aim points of an a, p. s. 
with spacing Zd between aim points. 

. Under assumptions (I) and (II) above, it follows that in Case 
] there are m hits if there is one hit and hence that 
P um, L) = P(n) ; that is, the value of P QU. m, L) 


is identical with the value of Р (п) for a simple salvo, For 


C, referring to Table 1, if 3 T/O's are aimed at al 
then three hits would be obtained with probabilities .250 and 
.383 , respectively, for the target position density functions 
f, and f, of Table I, Similarly, for a salvo of 6 T/O's, it is 
found that P,(3,2,L) has the values . 750 and „867 , respec- 
tively, for the target position density functions f, and f, of 
Table I, 

Turning to Case 2, the probability of more than one hit depends 


in part upon the amount of overlap!" between successive aim points 


of the a. p.s. To illustrate, consider a salvo of two T /O's fired 
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іп апа. р. з. ої two aim points with spacing 2d = 5 L. Thenfora 


target position density function £ xi) with ЈЕ = L (as in Table 1) 
1 


P (2, 1,5 Jasa? 597 


• 197 


Р.(2, 1,5 1) 

Finally, in Case 3, it is clear that if m T/O's are fired at 
each aim point of the a. p. s, with 2d < L , there is a possibility 
of m , 2m, etc. hits, Inthis case it is possible to calculate 
а 2d) where j is a positive integer having values l, 2, etc. 
up to some maximum value depending upon 2d, For a given salvo 
of nxm T/O's, it is important to note that whenever n> 1 and 
2d < L , a penalty is paid, namely, reduction in the value of 
Pun m, 2d) to achieve the probabilities Fan, m,2d) where 
nel, 

The spacing 2d = L has been measured along the x-axis, 
coincident with the assumed target course Со . Consider the 
T /O's on any track T, Other than normal to the x-axis. Let 
A denote the track angle between To and Со (measured as in 
Fig, 2 e), where for simplicity adjacent tracks are considered 
parallel). To intercept an interval L on the x-axis, the tracks 
must be spaced L sin A apart, L sin A is called the effective 
target length, For the one-dimensional case, it is clear that the 
probability of a hit associated with a given aim point does not 
depend upon the track angle. To avoid consideration of target 
width, very sharp track angles will not be considered here 


(roughly, A< 15° or A< 165° is not considered), 
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h(x, ;t) 


1/M 1/2M 
— — қ 
(и a 
— 
Figure 2 b: h(x, ;t) where Ex, - ML, М?1 
P(x,;t) 
~ ПИ 
E 1-N)L T 
(M+1+N)L 
Figure 2c: Р(х; ;t) where Ex, = ; Ex, 2 NL, M»N»l 
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n — — — — — —— Aim Point 
n=2 ل لا‎ Sequence 
nme س ي لل ل ر‎ (a. p. 8.) 
n= ——Ó—— SS 


Figure 2d: The a.p.s. depends upon n , the number of aim points, 


and the probability density associated with target position, 


L sin A 


T 
о 


Figure 2 e: Aim point spacing along the normal to T/O track. 
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CHAPTER 3 
ELEMENTARY PROBABILITY DENSITY FUNCTIONS FOR A 


STOCHASTIC VARIABLE WHICH DEPENDS LINEARLY ON TIME. ` 


Suppose that at a certain instant of time t , measurements 
are made of the position and speed of a point which is moving 
along a straight line at constant speed. Let the position and speed 
measurements both be subject to error and suppose that these mea- 
surements may be regarded as observations of two independent 
stochastic variables, Then it is frequently necessary to "predict" 
the position of the point at some later time t+ Y. In this con- 
text, "predict" means to determine the probability density func- 
tion associated with the position of the point at time t +T. 

The future position of the point is a stochastic variable which | 
depends linearly on time, If X denotes the future position of the 
point at time t + T and if x denotes the position of the point at” 
time t , then 

X2 x t sT 
where s is the speed, 

This chapter discusses the probability density function asso- 
ciated with X when the probability density functions of x and 
8 are (1) both uniform and (2) both normal. 

3,1 Probability Density Functions of x and s Both Uniform. 
Given: (a) A uniform probability density fi x { (xit) in 


I l 
the interval (x, > „Ex ‚x, + y Ex) associated with the 


> y 
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measurement of the position of the point on the line at time t . 
(b) A uniform probability density f, = £,(8;t) in the 

interval (в = iE, » 86 t EL) associated with the speed of 
the point on the line. When the position of the point at time t is 
known exactly and when time is considered as being measured 
without error, then the position of the point after an elapsed 
time TY is uncertain and the associated probability density 
function, here denoted by +) = fj(sT;t) , is uniform in the 
interval (s T - 3 ES Y, 827 + 5 Ез?) 

Clearly at Т = 0; f (xit) completely describes the pro- 
bability density associated with the point, with mean at хо 
and density fi z = in the interval (x, E j Ex, X + j Ex) 
and O elsewhere, This would be the case at any instant t 
when a measurement of x is made, lí no further measurement 
of x is available, but the position of the point at time 
t = t+ TY is required, eî it is necessary to consider the new 
stochastic variable 

Х = х+вт 

The probability density function f, of the new stochastic variable 


1 


is obtained by the convolution of £ ] and n . 


Thus, г = £3(X;t,) ғ.) Ж £,(X;t) 


3 


Eo - 87) f,(sT)rds 
"03S 1 
£X - x) f (x) dk 


го 
For the simple case here, the distribution diagram for f, is 
represented in Fig. 3a. 


l. Cf: M. E. Munroe, Theory of Probability, 1951, Section 38, 
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It is geometrically evident that: 
(1) The mean of X is moving, and at time t) isat x, t s T 
(2) The values of X are restricted to lie within the interval 


[x, + s T - (3Ex t J Es, x +s T + (¿Ex + 3 Es7)) 


which lengthens withT . 
(3) The area of the trapezoid representing the probability density 
function £, over this interval must always equal one (the total- 
ity of probability). Consequently, the ordinate values, which 
represent probability density, must necessarily be reduced as 


С increases, 


3,2 Probability Density Functions of x and s Both Normal, 


Consider now the case where each of the density functions f 


| 


f it^ is not uniform, but normal, In accordance with the nota- 


Z2 9 
tion in Appendix А 
== 5 
fi f (xx > 65%) 
f, git) 2 


== 


ҒЫТ” ар, fiit‏ = ر؟ 


uj 


SALILAT 
Then for any time с 93 Т, the stochastic variable 

X = х+вТ has a normal probability density function which, 

by the addition theorem for normal distributions may be repre- ` 


sented in the notation of Appendix A 


as f 


3 
where Хо x, + Bo T 


& 2065. fr 


x 


The consequent movement of the normal variate X , the 
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expansion of 0 X and the thinning out of f3 is also clear in this 


case as shown in Fig. 3 b. 
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£, (x, ) £ ($4 y 
ges M 2 | " 
1/Ex, 1/Es 
0 x 0 
хо ] о E 
f5(x) 


2a > ; 2b-2a — е 2b = Max[ Ex], Esr 
P ww «X _— 
х 


MUT 
Fig. 3a: f, = fi мб when fi and f, are 
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CHAPTER 4 
THE AREA OF UNCERTAINTY OF TARGET POSITION AND ITS 
ROLE IN THE SUBMARINE FIRING PHASE DECISION. 


4,1 Preliminary Considerations and Notation 
Let the independent stochastic variables (R,B) describe the 


position of a point ina plane; e.g., an observation of the range R 
and the true bearing B ofa target. Let the uncertainty of the 
exact value of R (at time t ) be associated with the uniform 

probability density function 

ГВ 1 1 

TER for R = JER SR Ss R, + 4 ER 

f (Rit) = 

0 elsewhere 


Let B be associated with the uniform probability density 


function 


1 1 1 
EB for В, - ЗЕВ< В < В, + 5 ЕВ 


£,(B:t) = 
0 elsewhere 


For definiteness, consider B measured in radians, say from a 
stationary submarine whose position will be considered a fixed 
datum point of the problem in any coordinate system desired, 
The joint probability density function associated with a pair of 
values (R,B) at t is then 

£¿(R, Bit) = f (Rit) (Ви) : 
Thus, f3 is positive over the annular sector shown in Fig, 4a 


as the domain D, = D,(t) 3 
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Suppose now that no matter what point in D, is selected, there 


3 
isa vector Y, uncertain as to direction and speed, that can send 
this point (R,B;t) into a new position, say (6, t, T) after an 

elapsed time ^ , As heretofore, consider T as stop-watch time 
measured from the instant t without error. As shown in Fig. 4b 


consider the angle 9 that \r makes with B to be an indepen- 


dent stochastic variable, with associated probability density 


elsewhere 


function 1 1 1 

Е5- for. 9. ~ z EÈ 5 0 < 0. + 5 ЕӨ 
Consider the speed S , i.e., the magnitude of АЛ, to be an in- 
dependent stochastic variable, with an associated probability 
density function 


1 
f- (S;t) = JES 


0 elsewhere 


1 Lx 
for 5, 7 27255 55 S, + „ES 


Then at the end of a specified time ^]" , the probability density 


associated with the magnitude of the vector — due to the variate 


S is, say b- l I. 
| EST for 36 = z ESY $ SY < 5. + > EST 
£.(STit) = | 
| 0 elsewhere 


Evidently, for each possible origin (R,B;t) in р, » the vector of 
uncertainty Y T has a domain of termination, say 
D, = р (в, В; АЛ, т st) consisting of the annular sector D, 


shown in Fig. 4b. And with D, is an associated probability density 


6 
function 


£,(S, 0; Tit) = £,(@st) £,(STst) 
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But V T may originate from any point in D, and further VT may 
assume any one of the set of values allowed by De . It may be consi- 
dered that ре "maps" D, into a new domain D- ж DAC A st T) 
where P and 8 are new stochastic variables, Functionally, D, 
may be represented as 
(4.1) О. = D, [R,£,(R):B,£,(B);8,£,(8);S,£,(S): Tt] 
The region D. may be obtained as the result of a simple geome- 
trical construction and is an eight-sided figure (see Fig. 4 c). 
The boundary of the area D. is comprised of the four sides of 
D; and the four sides of D, . With D, is associated a positive 
valued probability density function ғ. = t. (€, ЖТ)». Тһе 
mathematical problem is to determiine exactly how unit mass, 
representing the totality of probability, is spread over D; ғ. 
denotes the description of this mass distribution over р. 5 

No attempt to develop the analytic theory for the evalua- 
tion of the joint density function will be made here. Clearly, the 
probability density 1. of a set of values of the new variates | 
(0,6) at time (t+T) is zero at the boundary edges of D. and 
rises continuously, with possibly a different slope from each 
edge to a plateau of maximum probability density, say 

17 = max f, 

It is evident that there exists no simple way of describing the 
infinite variety of configurations of D, and the associated proba- 
bility density function f. for р. апа i depend upon the rela- 


tive magnitudes of all of the parameters indicated in formula (4. 1). 
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For example, Fig. 4 d shows the region D_ which results when D 


7 


the area of uncertainty of target position at time t , is propogated 


3 ? 


with no speed error and a discrete course error, 


As a point of reference in р. , 


of gravity of the mass distribution represented by Е. . The coordin- · 


define A, Yo be the center 


ates of A, will be designated (PP) or (0 o. uie T) when the 
time parameters are essential. When Ü and В ате independent 
stochastic variables then ч, = @ апа Во = f where e and В 
denote mean values of © and В respectively. 

Alternatively, if Е. has a maximum value at a unique point 
A, of р. ; then this point may be selected as a point of refer- 
ence of D, and in this case the coordinates of A, will be 
designated (О*, В") where 

f(P",B") = Max £06) . 
(6.8) € p, 


In the familiar special case where 1. is known (or approximated) 


in the form 


3 i (x-x ^ (у-у) 
EMT | О О 
, 7 ZT б; б, > : 

9 26 26 
^ у 


then A, = (х + Y): (See Appendix A, Fig. A-2). In this special 
case the iso-density contours defined by € =@ where € isa 
ее are concentric ellipses. 

Proceeding frcm the formulation of b. and f, with the 


simplest of assumptions concerning the component stochastic 


variates R,B,8,S, it would be desirable to turn to the case 
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where the component errors are independently and normally distri- 
buted. The reader will recognize the inherent mathematical diffi- 
culties which are beyond the scope of this paper. 
4.2 The Submarine Fire Control Problem 

The concepts and notation of Section 4. 1 will now be related to 
the submarine fire control problem, Thus, let (А,В) denote 
the actual range and bearing of MOT at the time of final observa-~ 
tion of the target (instant of last inputs into Position Keeper of 
TDC before firing). Ro and в. are the mean or expected yalues 
of R and B attime t with respect to the distributions f (Rit) 
and f,(Bit) , respectively, The functions f, and f, are 
assumed to be known, It will be assumed further that the values 
of R, and B o are identical with the final observations of R 
and B (inputs into TDC), Similarly, 90 and S denote the AOB 
(angle on the bow) and the speed of the target at ш» t. 8, and 
5. are the mean or expected values of 9 and S at time t with 
respect to the functions #4(9;+) апа #5(5;4) respectively. Again, 
the tons f 4 and f 5 are assumed to be known and the values 
of Ө апа So are assumed to be identical with the final obser- 
vations of € and S (inputs into TDC). Then the region D, 
may be called the area of uncertainty of MOT at the time of 
last observation oí the target. The region D. then represents: 
where the MOT may be found at time t+T (yunits of time after 
the last observation) when there are uncertainties in target 


range, bearing, course, and speed. The region D; will be 
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_ termed the area of uncertainty of MOT at t*'r. 

The interval ' may be considered as the sum of two intervals 

"Ті ала 12. Тат + 7, 

The interval Ti. represents the "dead-time" between the last 
observation of the target and the instant of firing; the interval 
To » which is a function of 7 js the time required for a "perfect" 
T/Oto reach a, = (0 B uit T) when the T/O is fired at the 
instant t+7,. Briefly, 7, is T/O running time and 
ET) 

The submarine is заа to have a TDC which can predict 
Ta within the interval т, i.e., T/O running time is calculated 
before the T/O must be fired. 

To this point, the direction of the target's motion has been 
discussed in terms of 9 , the AOB, It will be convenient to 
speak now of the target's course C , where, of course, 

C = C(@,B) and in particular C, - C(6, , B). Inasmuch as 
Со is the target's course as it appears in the TDC in accordance 
with previous assumptions, it is convenient to consider the aim 
points of the aim point sequence (see Ch, 2) of the T/O salvo to 
be along the course line Со extending through the point а. 

Unlike the simple one-dimensional сазе discussed in Chap- 
ter 2, there is in the two-dimensional problem no simple way 
to calculate the probability of a hit P: associated with the ith 
aim point of the a. p. s. when f. is generated by fi, 1, 5 


and f To discuss the two-dimensional case further, it is 


5 o 
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convenient to use the simplest of continyous probability density 
functions as a model for f Е 


To simplify them, consider 


1 Li 
ЕхҒу for х - Ех 5 х < 


1 
о о" pu 
Yo ZEY S y S y, t 38 


1A 


#2 = f(x,y) = | 
| 0 elsewhere 


With this simplification, DA? .B 5,7) js a rectangle which will be 
designated D(t+M). See Fig. 4 e. 

Recallthat D represents the area of uncertainty of position 
of MOT at time tt T. The time interval is computed for the 
collision instant of a "perfect" T/O with the point а’. Let C, , 
the "most likely" target course. through a, ‚ Бе identical with the 
x-axis of coordinates, Let the collision instant tty = to . То 
simplify further, assume that all T/O's have infinite speed through 
Dit) , and that all T/O's of a salvo pass through D(t,) at the 
same instant to ° The width of the target will be т белге ді For 


T/O's traversing paths parallel to the y-axis, the probability for 


a hit at a given aim point x, along C_ would be 





— 


х. + 7 L НИ ~ "A 
h(x,;t_) = f(x, y) dydx 
2' о 1 
=> = Ж L = бо 


which reduces to the result (2, 2) for the one-dimensional case of 
Chapter 2, For track angles other than 90? this is not the case. 
Referring to Fig. 4 f, it will be seen that the probability of a hit 
for such a track is simply the volume under f(x, y) "swept out" 


by the T/O "path" of width 1, зіп А . 
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Consider now the assumption that the T/O has infinite speed in 
the light of the simple example just discussed. Taking into account 
. the finite speed of the T/O, a very precise calculation of probabi- 
lity of a hit would require determination of the earliest possible 
instant, say (t, - b.) after firing at which the T/O enters a 
region D(t,- $1) in which f(t” бү) has a positive value and the 
late at possible instant, say (t+ 6,7 at which the T/O leaves a 
region D(t + 5.) in which f(t +62). has a positive value. When 
£(t, - 61) апа f(t +85) differ sufficiently over the "path swept 
out" by the T/O, there may be an appreciable error in assum- 
ing that the T/O passes with infinite speed through a region 
D(t,) and sweeps out a path under the probability density func- 
tion f(t.) : ÎÛ reaver, the statements above that the T/O 


, enters Es 6) and leaves D(t tó) imply that these regions 


have certain boundaries which will not be true when f has 
positive values over the entire plane as, for example, when at 
least one of the component probability density functions f L f, , 
f, » Or f; is normal, Under these circumstances it may be 
convenient to consider a region, say D(t,-5,) bounded by an 
€-valued iso-density contour of f(t - 5.) (cf, remarks follow- 
ing formula (3, 2)) which the T/O enters at t5791 and another 
region, say Dit,t 5.) , bounded by an €-valued iso-density 
contour of f(t + J which the T/O АЕ at ter 6, and deter- 
mine whether fit, - 5) and f(t +8.) differ "significantly" 


over the path swept out by the T/O. When these two functions 
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do not differ significantly, then the probability of a hit associated 
with an aim point may be visualized as the folume under КЕ; Е) 
swept out by the path of the T/O through D(t E) : 

In actual practice, the probability density functions fi, f; Р 
fa » and f; associated with target range, bearing, course, and 
speed inputs may not be known. However, it seems evident from 
the analysis of this paper that these functions are essential opera- 
tional data which are necessary for the submarine fire control 
problem if there is to be a rational basis for T/O salvos. It is 
likely that all four functions are normal density functions with 
standard deviations which may be estimated by operational trials. 

Appendix B gives a brief discussion of the numerical approxi- 


mation of f- by high speed digital computers when fi , f? : f, 
a 


and Е are given, 
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Figure 4 a: The annular sector D,(t) and the associated 
prabability density function f,(R, B;t) 
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Figure 4 b: The annular sector D A(R, B;v,T;:t) 
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Figure 4c: The region р. 
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Figure 4 e: The case where о. is a rectangle 


] in D 
f(D) =4 ExEy 7 
0 elsewhere 





Figure 4 f: The effect of T/O track angle on the area 


of D4 which is "swept out". 
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CHAPTER 5 
GENERAL REMARKS ON TORPEDO SALVOS 


бл T orpedo Lethality 


The.number of T/O's required, on the average, to destroy a 
target, depends on two primary factors: 

(1) The accuracy of the T/O, i.e. the probability the T/O 
will hit the target, and 

(2) The destructiveneas of the T/O, i.e. the probability 
that the target is destroyed if the T/O hits it. 

Here (1) is reduced to the problem of finding the Pis $ 
thence п. 

For (2) in many cases a single hit ig not enough to destroy a 
target. Even for merchant ships a single T/O hit is not usually 
enough to sink it, and heavy combatant ships are designed to 
withstand many T/O hits, To treat such cases exactly, one should 
determine the probabilities, say, di , d, ‘ d, , etc, oí sinking 
the ship if the ship is hit 1, 2, 3, etc. , times. Then if, for a 
given method of aiming T/O's, the probabilities of 1, 2, 3, etc. 
hits are P, : Р, i Р, » etc, the probability of sinking the ship 
would be di P, t d,P, + d,P, + ---, where the d's would repre- 
sent damage coefficients. 

For this paper it is simply assumed that there is a number m 
for each type ship such that, on the RARS m hits of that type will 


n 


eink the: hip with probability P, - Pím,n,d)- u P R 
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with d = 2 L sin A ;; the spacing along the normal to T/O tracks 


for the a. p. s. described in Chapter Z. 


5.2 Time Intervals in Salvo Fire 

Recall that, in firing torpedoes from a submarine, the firing 
is not a true spread in the sense that all T lo are launched at 
once to provide a simultaneous menace to the enemy. T/O's аге” 
fired singly and consecutively with a prescribed time delay to 
avoid countermining. | 

Let AY = time delay (an interval of deveral secs, ). Consi- 
der a spread of nxm T/O's, let gp- 1. Let T/ be initial 
make ready time after last inputs to TDC at time t. The mean 


n-e} 


T/O would go at time T,* AY after time t , onan 


expected run time =Y, say. Let mean Ti = is * ا‎ ay: 
then mean 77” =T ғт, and t = t+ 2 
The meantrack To , commenced at time t +7; may be 


considered as the datum line about which the spread, or aim 


point sequence (a. p. s.) is placed. 
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CHAPTER 6 


THE FIRING PHASE DECISION 


6.1 The Cardinal Rule 

The number of aim points n іп ап а, р. в. may be regarded 
as depending upon the function f(t ;€) introduced in Section 4, 2, 
Hence n is a function of У. | 

Therefore, asa cai dinal rule to minimize the spreading and 
thinning out of probability density, avoid protracted firings: 

FIRE AS SOON AS POSSIBLE AFTER DECISION TO 
COMMENCE FIRING PHASE. 

6.2 The General Procedure 

Most of the notation used here has already been discussed, 
otherwise it is defined in the glossary, Four sets of procedures 
are outlined, following the dictates of the four objectives pro- 
posed in Chapter 1. In each case n maximum i for which 
P. 2 pg and d = 7 L sin A along the normal to T/O tracks 
in the aim point sequence. Then P, = H Р. 

I, Fire Salvos to maximize — of at least one hit 
(i.e. case B1). Do not fire T/O's at aim points for which asso- 
ciated probability of hit is less than Po ° As nearly as possible, 
all T/O's are to be expended in not more than T days; obtain 
extra hits as indicated by frequency of attack. 


(a) # 2; € n fire one T/O at each of first $”, aim points 


(a. p. ) (cf. Section 2, 3) 


= Ші 
> Жү . : Inn ' POEM on 
} 
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(b) If n IH look ati then if: 
(i) 5] <%; fire J^ T /O's % each of first n a.p. and 
one T/O at each of first %; - An a.p. 
(11) n 52! < Du , fire حر‎ T/O's at each of first п а.р. 
(iii) XE < п < E fire one T/O at each of first n a.p. 
(c) When 2; = 22. request permission to reduce value of Ро 
(d) Return to base when all T/O's are expended, but do not 


remain on station more than T days. 


II, Fire T/O's to maximize probability of m_hits on targets 
requiring m hits to be sunk with probability Po (i.e., case Al). 


Do not fire at aim points for which associated probability of 
hit is less than Po ° 
As nearly as possible, all T/O's are to be expended in not 
more than T days; obtain extra hits as indicated by frequency 
of attack, 
(a) If mn 2), , fire m T/O!s at each of first \) а. р. 
(Optional to #те)}- mv T/O's at (Y+ 1)th a. p.). 
(b) I£ mn < %; » determine à; | then if: 
(i) м < у; fire J^. TfO's at each of first n a.p. 
(ii) mn< E < ў% fire А T/O's at each of first 
п а.р. 
(ii) F3: < mn < 5% fire m T/O's at each of 
first n a.p. 
(c) When 2j = у, request permission to increase 


value of Po and/or to reduce value of Po ° 
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(d) Return to base when all T/O's are expended, but do not 


remain on station more than T. days. 


III, Remain on station for ap proximately T, days, conducting 


as many attacks as possible. Ffre salvos to maximize probability 
of at least one hit (i.e., case B2) 


Do not fire T/O's at aim points for which associated probabi- 
lity is less than pg. As nearly as possible, expend all T/O's by 
end of patrol; obtain extra hits as indicated by frequency of attack, 

(a) If 5, $ n fire one T/O at each of first £i а. р. s 

(b Но < 221 fire Ач T/O's at each of first n a.p. 
and one T/O at each of first > - Ап а.р. . 

1 

(<) When 23; =z 2 Ж request permission to reduce value 
of pp. 

(d) Return.to base when all T/O's are expended, but make 


this as close as possible to Ti days. 


IV. Remain on station for approximately T, days, conduct- 
ing as many attacks as possible. Fire salvos to maximize probabi- 
lity of m hits on targets requiring m hits to sink with probability 


Poli ۰ case A2). 

Do not fire T/O's at aim points for which associated probabi- 
lity of hits is less than Po ° As nearly as possible expend all 
T/O's by the end of the patrol; obtain extra hits as indicated by 
frequency of attack. 

(a) If 251 5 mn fire m T/O's at each of first v, a.p. 
(Optional to fire 911 - mv, T/O's at (7) 1 +1) th a.p.). 
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(b) If mn < 25, fire AT /O's at each of first n a.p. (Optional 
to fire 231 - AyD T/O's at (n+ 1) th a.p.) | 

(c) When 24; = 257 request permission to increase value of 
Б and/or reduce value of Ро. 

(d) Return to base when all T/O's are expended, but make this 


as close as possible to Ti days. 
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Appendix A 


The Normal Distribution 


A-1 The One-Dimensional case. 

The normal distribution is one of the most significant mathema- 
tical models, from both the theoretical and practical standpoints, 
upon which the techniques of applied statistics are founded. ' 

Using as an example, say, a set of n measurements, with 
errors x, (i ='1,2,--2n), some underlying assumptions of the 
distribution are; 

(a) In. performing the set of measurements, not all the errors 
x, are equally likely to occur. 

(b) In general, large errors are less likely to occur than 
small ones, 

(c) If there is a best value m , then positive and negative 
errors about this value are equally likely to occur, 

(d) If (c) is true, then m is the arithmetic average of the 
individual measurements, 


The density of this distribution is given by 


tag BG) 


Thus, a normal distribution is completely determined by specify- 





ing its mean m , and standard deviation Ó i.e. f(x:m,( ). 
The function is shown in Fig. A-I for several values of 6 . 
Changing m merely shifts the curves to the right or left without 


changing their shape. 





By the linear transformation y = — , Қх) is put in stan- 


dard form, with m = ovand 6 = 1. The density function is then 


{х;0,1) = ap Exp ( - x ) 


It is in this form a. usually finds the integral tabulated for 
any value of x. 
A-Z The uou cnsional'CA&E. 

A simple case of two independent, normally distributed varia- 
bles will be discussed here, Given fi(xix,, 6:2) and EYY? 2 
such that f(x,y) = f, (x) f,(y) . Let the point (x. yg) = а. 


The joint density is given by: 2 2 | 
1 1 imi. %-%5 
2 2 


(i.e. unit mass is spread over a symmetrical mound). Fig. A-2, 





Curves formed by a vertical slice parallel to an axis are of 
the same shape as the normal distribution curve. 


The curve formed by a horizontal slice is an ellipse 


2 2. 
х-а у-а 
f(x,y) = Ө + (ғ e) қ = (a constant) 


The equation represents an ellipse with center at a, semi- 








axes of length (c б, с ©) and axis parallel to (x,y) axes. Ап 
ellipse of this nature is called an iso-density contour ellipse. 

The probability that (x,y) lies inside the ellipseequats P when 

E =X 4 » f = 2 (where = is found using the x,y values and 
P is the probability in per cent for the corresponding fractile of 


the Chi-squared distribution with 2 degrees of freedom). 





To define a bound for the area D covered by f(x,y), any 


contour ellipse € may be selected, e.g. noting that: 


2 
NZ 
- 4.61 (f=2) 
А 


then ae = 4.61 8 = 2.15 


let € = 90% contour ellipse with semi-axes (2. 156. , 2; Po) 


i.e. for E = 90% the ellipse is 4.30, by 4. 36 د‎ centered at a. 





Figure A-1: Normal Frequency Curves. 
m=0; 6 = 0,4,1.0,2.5 
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Figure A - 2: Bivariate Normal Distribution 
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Appendix B 
Proposal for Preparation of Histograms Used in 


Salvo Firing 


A proposal for using an electronic digital computer for deter- 
mining D and F(D) et for use as reference diagrams, 
First, prepare statistically, the histograms representing 
R, (В); В, (В); Ө, f(9) where 9 = 0(C,B); and S, f(S). 
Thus for R, we might have the histogram shown in the following 


figure. 





For each component the machine is instructed to pick a random 
number (by a shb-routine method) say 0 < r < 1 and this num- 
ber instructs the machine from which step oí the cumulative graph 
to pick the value of R. e.g.,0 S r & A, picks R, А, < rs 
A, " A, picks R5, etc. In this manner, a value of each compon- 
ent is picked and the machine computes the simple trigonometry 
involved to get the X, Y coordinates, say from the datum point 


for a given т. 


m 





These are compared to the grid values, sorted, and added to the 


proper cell as a simple count of.the occurrence 


= x3 х3 4 х; 
Уі Г ЕЕ 5 
„| | o [t2 
I OI ж: 
3 100 yard 
1 
жиектер БЕРЕ 


D and f(D) for various set of components and J's could 
be make up to guide the approach officer in the selection of the 
р, 5 and n . 

Not only would one like histograms based on various mean 
values to use in the computation, but also on the method of 
measurement, such as range histograms at various ranges 
taken by: radar, optical, sonar, etc: Bearing: optical, hydro- 
phone, sonar; Speed; screw count; tracking, etc. 

lf each histogram was made up in standard steps the compu- 
ting routine could be used over and over again to construct D's 
and f(D) for various combinations by just changing the values 
in the cumulative histogram cells. 

Note that the computer would produce discrete mounds that 
would be as good as the histogram data from which they were 


prepared. 
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